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Abstract. Physical constraints such as positivity endow the set of quantum 
states with a rich geometry if the system dimension is greater than two. To 
shed some light on the complicated structure of the set of quantum states, we 
consider a stratification with strata given by unitary orbit manifolds, which can 
be identified with flag manifolds. The results are applied to study the geometry 
of the coherence vector for n-level quantum systems. It is shown that the unitary 
orbits can be naturally identified with spheres in R" ~^ only for n = 2. In higher 
dimensions the coherence vector only defines a non-surjective embedding into a 
closed ball. A detailed analysis of the three-level case is presented. Finally, a 
refined stratification in terms of symplectic orbits is considered. 
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1. Introduction 

The Bloch vector |jy provides a representation of the quantum states of a two-level 
system in terms of real observables, and allows the identification of quantum states 
with points in a closed ball in 3D Euclidean space, the Bloch ball, which has proven 
to be extremely useful. In quantum information theory, for instance, the states of 
a single qubit can be identified with points on the surface of the Bloch ball (if the 
state is pure) or points inside the ball (if the state is mixed). Unitary operations can 
be interpreted as rotations of this ball, and dissipative processes as linear or affine 
contractions of the Bloch ball ^ |2| . 

Naturally, many efforts have been made to generalize the Bloch ball to higher 
dimensions by defining a generalized coherence vector similar to the Bloch vector for 
two-level systems. However, while it is easy to define a Bloch or general coherence 
vector ^ Ej for A'^-level systems, it has become clear that the geometry of quantum 
states in higher dimensions is far more complex than in the two-level case. Some efforts 
have been made at determining the set of Bloch vectors corresponding to physical 
states for higher dimensional systems using the higher trace invariants or Casimir 
invariants, for instance by Kimura 6 , who also demonstrated the complicated and 
asymmetric structure of the set of Bloch vectors in higher dimensions. Byrd and 
Khaneja "7^ independently provided a similar characterization of the positivity of the 
density matrix in terms of the coherence vector representation. 

In this paper we pursue a different approach to study the structure of quantum 
states for higher dimensional systems, and the origin of the difference between the 
two-level case and higher dimensions. In section II we define a natural stratification 
of the set of density matrices in terms of unitary orbits. We show that the unitary 
orbits can be identified with flag manifolds and determine their dimensions. For two- 
level systems this stratified set can be identified with a smooth real manifold with 
boundary, the 3D Bloch ball, with strata given by concentric spheres. 

In section III we briefly review the definition of the Bloch vector and define a 
general coherence vector, which can be naturally embedded into a closed ball in E" ~^ . 
This embedding is surjective for two-level systems, hence allowing the identification 
of physical states with points in this closed ball. We show that for higher dimensional 
systems, however, the embedding is never surjective due to positivity constraints, 
and the dimensions of the orbit manifolds vary depending on the multiplicities of the 
eigenvalues of the states in each stratum. Furthermore, for two-level systems there 
exists a total ordering of the strata given by the length of the Bloch vector, or the 
distance of the unitary orbit from the center of the Bloch ball ~ the Bloch vector of 
pure states has length one, and the shorter the Bloch vector, the more mixed the 
state is. For higher dimensional systems we show that the length of the coherence 
or generalized Bloch vector is no longer a sufficient measure for the disorder of the 
system. A detailed analysis of the three-level case concludes section III. 

Finally, in section IV we briefly consider a refined stratification of the set of 
density matrices defined by the action of the symplectic group. The symplectic group 
is of interest because it occurs naturally for quantum systems with certain dynamical 
symmetries (such as atomic systems with degenerate energy levels) and it is the only 
proper subgroup of the unitary group that acts transitively on certain unitary orbits 
[TO] . We show that the symplectic orbits of pseudo-pure states agree with the unitary 
orbits and provide bounds on the dimensions of other symplectic orbits, showing that 
the symplectic orbits generally have much lower dimension. From a control point of 
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view this means that a 2n-dimensional Hamiltonian control system with dynamical 
Lie group Sp(n) is pure-state controllable but we cannot control generic ensembles. 

2. Density matrices and unitary orbits 

Throughout this paper we restrict our attention to quantum systems whose Hilbert 
space is a finite-dimensional complex vector space C", where n is an integer greater 
than one, and for simplicity we will often use V to denote C" with the standard 
Hermitian inner product (•,•). Any physical state of this system can be represented 
by a density operator, i.e., a positive semi-definite (self-adjoint) linear operator with 
trace one. The subset of rank-one density operators corresponds to pure states of the 
system, all other density operators to mixed states. 

In the following we denote the set of all positive semi-definite operators on V by 
P(V^), the subset of all density matrices by D{V), and the subset of pure states by 
Diiy). We also define the class of pseudo-pure states or pure-state- like ensembles as 
the set of density operators Dy{V) whose spectrum consists of exactly two distinct 
eigenvalues with multiplicities one and n — 1, respectively. It is well known that D{V) 
forms a convex subset of the set of Hermitian matrices since given two density matrices 
pi,P2 € D{V), the straight hue path defined by T{t) := (1 — t)pi + tp2 for t € [0, 1] 
is contained in D(V). To see this, note that T{t) is Hermitian, Tr(r(t)) = 1 and 
{T{t)x \x) = {1- t){pix I x) + t{p2X \x)>0 for all a: G C" and t G [0, 1]. 

We define the orbit of a quantum state p under the action of the dynamical Lie 
group G to be the set G • p := {g ■ p ■ g~^ \ g G G}. The orbits endow the set of 
quantum states with the structure of a stratified set. In principle, any Lie group G 
acting on the set of density operators defines a stratification of D{V). However, the 
orbits under the action of G = U(n) are of particular interest since they determine 
the most general evolution of the quantum states in a closed system. From a control 
point of view, the unitary orbit of a state represents the maximal set of states that are 
reachable from the given state via (open-loop) coherent control, or in the language of 
quantum computing, by applying a unitary gate to the state. 

Before we attempt to classify the orbits, we recall the following standard result 
from linear algebra: 

Proposition 1 Let pi and p2 be two density matrices. The following are equivalent: 

(i) pi and p2 are unitarily equivalent, i.e., p2 = UpiU^ for some unitary matrix U. 

(a) pi and p2 have the same spectrum, i.e., the same eigenvalues including 
multiplicity. 

(Hi) Tr{pl) = Tr(p5) for all r = 1,2, ... ,n. 

This result shows immediately that the orbit of a density matrix under U(n) 
is uniquely determined by its spectrum, i.e., two density matrices belong to the 
same unitary orbit if and only if they have the same eigenvalues with the same 
multiplicities m. Each orbit can therefore be represented by a canonical diagonal 
density matrix of the form 

p = diag(Ai/„i, . . . , ArI„J, (1) 

where the eigenvalues shall be ordered such that Aj > Xj for i < j to ensure a unique 
representation. Since the Aj can be arbitrary real numbers in [0, 1] provided they 
satisfy X^^^i Aj = 1, we see immediately that there are infinitely many distinct 
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orbits corresponding to the (uncountably) infinitely many possible choices for the A^. 
Thus, we can say that the unitary group U(n) partitions the set of density matrices 
D(V) into an uncountably infinite family of (distinct) orbits or strata. 

We can define a (partial) ordering on this stratification via majorization. Let pi, 
P2 be two density operators with eigenvalues Om , i = 1,2, counted with multiplicity 
and ordered in a nonincreasing sequence, pi -< p2 if 

X]m=l — Sm=l ^rn , /c = 1, . . . , 71 — 1 , , 

^(1) _ (2) 
Z^m=l "m — Z^m=l "™ ' 

For instance, pi = ■|diag(l, 1, 3) ^ p2 = ■5diag(2, 2, 1) since 1<2, l + l<2 + 2 and 
1 + 1 + 3 = 2 + 2 + 1. Majorization has been shown to be a useful way to compare 
the degree of disorder of physical systems W and naturally defines a partial ordering 
on the unitary orbits (strata) via 0i ~< 02 if 0i — 0[pi] and pi ^ p2- However, note 
that only some orbits can be compared that way. Consider pi — idiag(5, 2, 1) and 
p2 = |diag(4, 4, 0). We have 5 > 4 but 5 + 2 < 4 + 4. Hence, neither orbit majorizes 
the other. 

To determine the nature of the strata given by the orbits, we define the isotropy 
subgroup or stabilizer at p as the subgroup Gp of elements in G that leave p invariant, 
i.e., for which we have g ■ p ■ = p. We shall show that the orbit of an element 
p G D{V) under the unitary group U(n) can be identified with a certain type of 
manifold called a flag manifold. For the purpose of proving this result we observe that 
U(n) is a compact Lie group and hence a compact topological group| and the space 
of density matrices is a Hausdorff space§, and we have the following result (see [S], for 
instance) : 

Proposition 2 If G is a compact topological group acting on a Hausdorff space X 
and Gx is the isotropy group at x then the map (j) : G/Gx ^ G-x is a homeomorphism. 

Theorem 1 Let U(n) act on D{V) by conjugation and let p be a quantum state with 
r > 1 distinct eigenvalues Xi with ( geometric ) multiplicity n,; . Then the orbit of p is 
homeomorphic to the flag manifold 

U(n)/[U(ni) X U(n2) x • • • x U(n^)] 

of real dimension — X]i=i 

Proof: Let Ei be the eigenspaces of p with dimi?.; — Ui. Since p is unitarily 
equivalent to the diagonal matrix diag(Ai/„j , . . . , Xrln^), we have an orthogonal direct 
sum decomposition ofV — C" of the form V — Ei (B ■ ■ ■ (B E^- g & U (n) stabilizes p if 
and only if g preserves the eigenspaces Ei, i.e., the restriction of g to each eigenspace 
must be an isometry, i.e., g preserves the eigenspaces Ei. Hence, the orbit of p is 
homeomorphic to the flag manifold U(n)/[U(rti) x U(n2) x • • • x U(nr)] by Proposition 
m ■ 

I A topological group is a topological space X endowed with a group structure that allows us to 
"multiply" elements of the space and compute inverses such that both operations are continuous with 
respect to the topology. The unitary group, for instance, is a multiplicative group since multiplication 
of two unitary matrices gives a unitary matrix, every unitary matrix has an inverse given by the 
Hermitian conjugate, and the identity provides a neutral element. Furthermore, as a subset of 
the complex matrices the unitary group is naturally endowed with a topology that allows us to 
separate two unitary matrices by open sets, and matrix multiplication and Hermitian conjugation 
are continuous with respect to this topology. 

§ A Hausdorff space basically is a space endowed with a topology that allows us to separate points 
by disjoint open sets. 
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Corollary 1 // p G Dii{V) (pseudo-pure state) then the orbit of p is homeomorphic 
to U(n)/[U(1) X U(n — 1)], which is homeomorphic to the complex projective space 

To illustrate the result, we explicitly compute the orbits under the action of U(n) 
for n = 2 and n = 3. 

Example 1: Let p be a 2 x 2 density matrix, p is unitarily equivalent to 
diag(r, 1 — r) with 0<r<l. Ifr=l — r then p — i/2 and the orbit of p is 
homeomorphic to U(2)/U(2), i.e., a single point. Otherwise, its orbit is homeomorphic 
to U(2)/[U(1) X U(l)] ~ CP^ 

Since CP^ is diffeomorphic to the sphere S^, this shows that any U(2) orbit of 
a two- level system is homeomorphic to S*^, except the trivial orbit of the completely 
random ensemble which consists of a single point. 

Furthermore, note that the requirement of positivity of p reduces to < ^ < 1 and 
hence < r, or equivalently, Tr(p^) = + (1 — r)^ = 1 — 2r + 2r^ < 1 — 2r + 2r = 1. 
We shall see that this implies that the set of all 2 x 2 density matrices (the union of 
all orbits) is homeomorphic to a closed ball in R^. 

Thus, we have a neat mathematical justification for the Bloch ball description of 
a two-level system, which will be discussed in detail later. 

Example 2: Let p be a 3 x 3 density matrix. If p has only one eigenvalue with 
multiplicity 3 then p = ^I^ and its orbit is homeomorphic to U(3)/U(3), which is a 
single point as before. 

If p has two distinct eigenvalues then it is a pseudo-pure state unitarily equivalent 
to p = diag(l — 2a, a, a) where < a < 1 and a 7^ i. Its isotropy subgroup is therefore 
U(l) X U(2), and its orbit is homeomorphic to U(3)/[U(1) x U(2)] and has dimension 
9 - 1 - 4 ==4. 

If p is a generic ensemble with three distinct eigenvalues a, b, c then its canonical 
form is p = diag(a, 6, c) and its isotropy subgroup is U(l) x U(l) x U(l). Hence, its 
orbit is homeomorphic to U(3)/[U(1) x U(l) x U(l)] and has real dimension 9 — 3 = 6. 

The results of the previous example are summarized in the Tabled The table 
also provides a complete classification of the orbits for n ~ 4. 

The previous two examples clearly show the difference between two-level systems 
and higher dimensional systems (n > 2). While all orbit manifolds (except the trivial 
orbit of the completely random ensemble) for two-level systems are homemorphic 
to a sphere, no such homeomorphism is possible in the latter case, i.e., the orbit 
manifolds for higher dimensional systems can never be identified with spheres in a 
higher dimensional Euclidean space. That is, although we can always embed the 
quantum states of the system in a compact subset (closed ball) of a real vector space of 
sufficiently high dimension, there is no one-to-one correspondence between spheres in 
this ball and orbits of quantum states, except for n = 2, which explains the difhculties 
one encounters when trying to generalize intuitive reasoning valid for the Bloch ball 
for n = 2 to Bloch vectors in higher dimensions. 
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n = 3 


manifold 


dim. 


p = diag(a, a, a) 
p = diag(a, b, b) 
p = diag(a, b, c) 


point 

U{3)/[S' X U{2)] 
U{3)/[S^ x 51 X 5*1] 




4 
6 



n = 4 


manifold 


dim. 


p = diag(a, a, a, a) 
p = diag(a, b, b, b) 
p = diag(a, a, b, b) 
p = diag(a, b, c, c) 
p = diag(a, b, a, d) 


point 

UiA)/[S' X C/(3)] 
U{A)/[U{2) X U(2)] 
U{A)/[S^ X 51 X U(2)] 
U{4-)/[S^ X 51 X 51 X S^] 



6 

8 

10 
12 



Table 1. Manifolds and their dimension for the unitary orbits of quantum states 
based on their canonical form. All parameters a,b,c, . . . in the table above are in 
[0, 1] such that Tr(p) = 1 and different letters represent different values. 



3. Coherence vector and embeddings of quantum states 

3.1. Definition of Block or cokerence vector 

For a two-level system any density operator can be expanded as 

P= ^(-^2 + Xaj; + y<7y + ZCTz) 

where h is the 2D identity matrix, and = |1)(2| - |2)(1|, Oy = i(|l)(2| + |2)(1|) and 
(Jz = i(|l)(l| — 12)(2|) are the usual (unnormalized) 2D Pauli matrices. The coordinates 
X, y and z are real since the Pauli matrices are Hermitian and * = Tr(/9cr,) for 
* = x,y, z. Hence, the state of any two-level system can be characterized completely 
by the real vector s = (x, y, z), called the Bloch vector. 

For an n-level system we can proceed in a similar fashion. Let 

<,s = \r){s\ - \s){r\, al, = i{\r){s\ + \s){r\), 

ol = i^^ (ELi l^> - + 1) + 1 1) 

for 1 < r < ri — 1 and r < s < n be the generalized Pauli matrices in dimension n. 

2_-| 

The set \pkW^\ = {c^r s> crj^s, 1 1 ^ ^ < ''' < s ^ 't-} forms a basis for the space 
of n X n traceless Hermitian matrices satisfying the ortkogonality condition 

(c^fe I ^i) = Tr(CTfeCT^) = 2Sk,e- (4) 
Every density matrix can be expanded with respect to this basis 

1 1 "'"^ 

P = -In + o X] ^'^^'^ 
71/ ^ 

fe=i 

where Sk ~ Tr{pd-k) for fc = 1, ... — 1. The resulting real vector s = {sk)^^^ is 
the Bloch vector of the n-dimensional system. 

Although the Bloch vector defined above is useful, it is generally more elegant, and 
often more convenient, to work with an orthonormal basis. To this end, we define the 
normalized Pauli matrices {ak}^J[^ = {:^^fc}fe=i^' which satisfy the orthonormality 
condition 

(o-fe I <Jt) = Tv{ak<Jt) = 5k,t. (6) 



Orbits of quantum states and geometry of Block vectors for N -level systems 



7 



Furthermore, {a^ \ k — 1, . . . , n'^ — 1} together with ctq ~ "^-^n forms an orthonormal 
basis for all Hermitian nx n matrices, and we can expand any density matrix in terms 
of this ON basis 



where Sk = Ti{pak) for A; = 0, . . . , — 1. Since 1 — Tr(p) ~ \/nso, we have sq = 

for all density matrices. Hence, p is completely determined by the real — 1 vector 
s = (si, . . . , Sn2_i). This vector is often called the general coherence vector. Equations 
(|SJ| and O are equivalent, and we easily see that s = 2s, i.e., the standard Bloch vector 
differs from the coherence vector only by a factor of 2. 

It is easy to verify that the mapping that takes p to the real coherence vector s, 
or equivalently the Bloch vector 5 = 2s, defines an embedding of the density matrices 
into a closed ball in M" ~^ for all n > 1. However, the two- level case is special in that 
the embedding defined is surjective. We shall now discuss the nature of the resulting 
differences between the n = 2 and n > 2 case, and provide a detailed analysis of the 
three-level case. 

3. 2. Bloch ball picture for n = 2 

For two-level systems the embedding defined above is not only one-to-one but also 
surjective, and hence provides a homeomorphism between orbits of density matrices 
under U(n) and the closed ball of radius one (Bloch vector as defined above) or radius 
i (coherence vector as defined above) in M'^. Unitary transformations of a density 
matrix can be interpreted as real rotations of the this ball. Example ^ shows that 
the action of the unitary group on the set of quantum states partitions it into an 
uncountably infinite number of distinct orbit manifolds, homeomorphic to concentric, 
two-dimensional spheres in the Bloch ball, with the exception of the trivial orbit of 
the completely random ensemble, which is mapped to the single point at the center of 
the ball. 

It is also easy to see that the distance of an orbit from the center of the ball is 
determined by Tr(/9^) via f = 2Tr(p^) — 1 (Bloch vector) or r ~ Tr(/9^) — ^ (coherence 
vector). Pure states (Tr(p^) ~ 1) have maximal distance from the center of the ball 
and hence form its boundary. Furthermore, the disorder of an orbit is completely 
determined by its distance r(0) from the center, i.e., 0i 02 if r{0i) < r(02). 

From the point of view of controllability of quantum systems, it is also worth 
noting that all orbits (with the exception of the completely random ensemble) have 
the same dimension and geometry. Hence, any group that acts transitively on the 
class of pure states, for instance, will also act transitively on all classes of mixed states 
and vice versa. Of course, the only such groups are U(2) or SU(2). Hence, pure-state 
and mixed-state controllability are equivalent notions for two-level systems. 

3. 3. Bloch ball picture for n > 2 

For n > 2 the situation is quite different due to the fact that the embedding into 
a closed ball in K" ~^ defined by p i-^ s is not surjective. It is easy to see that the 
distance of each unitary orbit from the center in M" remains completely determined 




(7) 



k=0 
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by Tr(p2): 

1 

iNr- E «^ = TV)-i- (8) 

k=l 

However, a glance at the dimensions of the orbits in Tabled shows immediately that 
each orbit is only a submanifold of a sphere of a fixed distance from the origin. For 
instance, as we have shown in the previous section, the orbit of pure states for n = 3 
corresponds to a four-dimensional submanifold of the seven-dimensional boundary 

sphere with radius ^1 — ^ in R^. Since there is only a single orbit of pure states, the 
remainder of the points on the boundary sphere do not correspond to physical states. 

For mixed states the situation is more complicated since each sphere of fixed 
radius r from the origin now generally contains an uncountably infinite number of 
distinct orbits, all satisfying Tr{p„i) = 1 and Tr(p^J = but differing in higher 
trace invariants, e.g., Tr(p^) ^ Ti{p^,) for some 2 < k < n. The dimensions of the 
orbits contained within each sphere vary depending on the type of ensemble but each 
sphere (except for the boundary) generally contains a set of positive measure of orbits 
corresponding to physical states, and often a positive-measure set of points that do not 
belong to physical orbits. Moreover, the degree of disorder of an orbit can no longer 
be properly characterized by its distance from the center of the ball. Orbits contained 
in the same sphere can often not be compared with respect to our partial ordering and 
may have different von-Neumann entropy. Furthermore, orbits at different distances 
from the origin may have the same von-Neumann entropy. 



3.4- Analysis of n — 3 case 

To illustrate the general statements above, let us consider the ri = 3 case and the set 
of orbits determined by the family of states p = diag(a, 6, c) with b, c given by 

b + 

c =^{l-a-K), 



K = \/-l + 2a-3a2 + 2c2. 

Note that we have Tr(/9) — a + b + c — 1 and Tr(/9^) = C2 for all a. However, for 
p to represent a physical state a,b,c must be real and have values in [0,1]. These 
constraints imply that the argument of K must be non-negative and K < 1 — a, or 
equivalently < (1 — a)^. This yields the inequalities 

3a^ - 2a+ 1 < 2c2, (9) 

2a^ -2a+l>C2, (10) 

which must be simultaneously satisfied. To ensure that there is a one-to-one 
correspondence between parameter values (0,02) and orbits, we further require a > 
b > c. Since a + b + c = 1 is implies a > i. The constraint & > c is automatically 
satisfied because K is real and > 0, whereas the constraint a > b implies: 

6a2-4a+l>C2. (11) 

Inequality ((HJ is satisfied for 

l + Ki' 

a e 
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Figure 1. Diagram indicating parameter values corresponding to piiysical orbits 



for Ki = V6c2 — 2. Inequality is satisfied for 



a e 





u 


■f - 








3 'J 



For C2 < ■k inequality (|10|l is satisfied for a G [0, 1], 



and C2 > i it is satisfied for 



a e 



0, 



I + K2 



for K2 = \/2c2 — 1. Combining these inequalities and noting that a > 



I leads to 



C2 < 



C2 > 



a e 



a e 



^1/6), J(l 



^1) 



(12) 
(13) 



2a 4 
- 2a 



2C2, 
= C2. 



Fig. n illustrates the situation. The solid curve corresponds to 3a^ - 
the dash-dot line to 6a^ — Aa + 1 — C2 and the dashed line to 2a^ 
The points {a, 02) on the sohd hne correspond to orbits of pseudo-pure states p — 
diag(a, 5(1 — a), 5(1 — a)), and the points on dash-dot line correspond to orbits of 
pseudo-pure states of the form p — diag(a, a,l — 2a). The points below the solid 
curve correspond to non-Hermitian matrices. The points (a, C2) above the dashed 
line correspond to non-positive Hermitian matrices. The points between these two 
curves represent physical states. However, only points in the shaded region between 
the curves satisfy all inequalities and represent unique orbits. 

The figure also shows that for any C2 £ (51 1); there exists a positive-measure 
set (interval) of a-values that correspond to distinct physical orbits, embedded in a 

sphere of fixed radius r = \J C2 — \ from the origin. This shows that the number of 

distinct orbits contained within each sphere of radius < r < is uncountably 
infinite. A unique orbit of fixed distance from the center of the ball exists only for 
the special cases C2 = ^ and C2 = 1, the former corresponding to the trivial (i.e., 
zero-dimensional) orbit of the completely random ensemble a = 6 = c = ^ that forms 
the center of the ball, and the latter to the four-dimensional orbit of pure states 
contained within the boundary sphere of the ball. Most of these orbits are generic 
and hence have dimension six. However, each sphere also contains at least one (two if 
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Figure 2. a + 6 = (l + ci + K)/2 as a function of a for various values of C2 



C2 < 5) pseudo-pure orbits of dimension four, corresponding to the points (0,02) on 
the boundary curves 2c2 = 3a^ — 2 + 1 and C2 — 6a^ — 4+1, respectively. 

Furthermore, for C2 < 5 all values of (a, C2) that correspond to Hermitian matrices 
actually correspond to positive Hermitian operators, i.e., physical states. Hence, the 
union of all orbits with C2 < ^ fills a ball of radius For C2 > ^, however, the 
positivity constraint kicks in and eliminates more and more a-values as C2 approaches 
one. This means that the physical orbits inside the ball get sparser as C2 increases. 
However, for each C2 < 1 the orbits always occupy a positive measure set of the sphere 
they are embedded in since for each C2 < 1 there is a ID set of positive measure 
of a values representing disjoint physical orbits of dimension six. Hence, the union of 
these orbits (not counting the one or two pseudo-pure orbits) occupies a set of positive 
measure inside each seven-dimensional sphere 5^. For C2 = 1 the positivity constraint 
eliminates all a-values but a — 1 and hence the boundary sphere contains only the 
single four-dimensional orbit of proper pure states. 

As regards the ordering of orbits, we have 0[a,C2] ^ 0[a, Cj] if C2 < c'2 since 
K{a,C2) < if (0,03) for C2 < C2 and a fixed. In general, however, we cannot compare 
orbits that have the same distance from the origin since a + b = ^(l + a + K) is 
monotonically decreasing over the valid range of a for C2 > ^ , and non-nionotonic for 
c < i, i.e., a < a' typically implies a + b > a' + b'. See Fig. [21 for a plot of a + 6 = 
^{l + a + K) as a function of a for various values of C2. Fig.O for comparison, shows a 
plot of the von-Neumann entropy of the orbits 5 [0] = — [alog(a) + 61og(6) + clog(c)] as 
a function of the largest eigenvalue a for various values of C2 . Note that for C2 > ^ fixed, 
the von-Neumann entropy increases as function of a, for c < ^ it is non-monotonic as 
function of a. Also note that for sufficiently large values of S, there are many orbits of 
varying distance from the center with the same von-Neumann entropy S although for 
fixed a the orbits with the largest distance from the center have the smallest entropy, 
which is in accord with our partial ordering of the orbits. 

4. Actions of the symplectic group 

In the previous sections we have shown that the action of the unitary group on the 
set of quantum states endows it with the structure of a stratified set and studied the 
properties of the strata defined by the unitary orbit manifolds. This stratification of 
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Figure 3. von-Neumann entropy of orbits as a function of a for various values 
of C2 



the set of density matrices was justified on physical grounds since the dynamics of a 
(closed) quantum system is usually determined by the action of the unitary group. 

Some physical systems, however, exhibit symmetries that restrict the dynamics 
to a proper subgroup of the unitary group. For instance, the dynamics of an atomic 
system comprising two n-fold degenerate energy levels subject to coherent control fields 
of various polarizations is governed by the symplectic group Sp(n) due to dynamical 
symmetries |1()| . Physical systems that exhibit symplectic symmetry are also of special 
interest in quantum control since they are pure-state controllable but lack mixed-state 
controllability JU]. This means, for instance, that we can steer the system from any 
pure initial state to any other pure state; however, if the system is initially in a mixed 
state, then it may not be possible to steer it to another mixed state even if this state 
is unitarily equivalent to the initial state, since the set of unitary operations at our 
disposal is limited and there may not be a sympletic unitary transformation that 
achieves the desired aim |11) . 

The action of the symplectic group also induces a stratification of the set of density 
matrices. Since the sympletic group Sp(n) is a subgroup of the unitary group U(2n), it 
naturally follows that the sympletic orbits are smaller than the unitary orbits. Hence, 
there will be more sympletic orbits. We can also think of the action of the symplectic 
group as partitioning the unitary orbits into symplectic suborbits. The stratification 
induced by the symplectic group is therefore a refinement of the stratification induced 
by U(2n). The remainder of this section is devoted to exploring the relation between 
the symplectic and unitary orbits / stratification. 

Mathematically, the symplectic group Sp(n) is the subgroup of unitary 
transformations A E \J{2n) that satisfy A^JA = J for 



where /„ is the identity matrix in dimension n. Note that, technically, any group G 
that satisfies {S'^ JS — J \ V5 G G} for a matrix J unitarily equivalent to the J in 
H14|l is a representation of Sp(ri) but we shall assume the standard representation with 
J as in H14|l here. The condition S'^JS = J then implies that any 2n x 2n complex 




(14) 
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matrix S G Sp(n) must be of the form 

S=[_% jl) (15) 

where A and B are n x n complex matrices and A* denotes the complex conjugate of 
the A. 

We now show that Sp(7i) acts transitively on the unitary orbits of some states, 
but for the majority of states the symplectic orbits have lower dimension than the 
unitary orbits. The following results are an extension of earlier results showing that 
Sp(n) acts transitively only on pure-state-like and completely random ensembles |11| . 

Proposition 3 If p is the completely random ensemble ^l2n or a pseudo-pure state 
then its orbit under Sp(n) is the same as the orbit under U(2n), i.e., Sp(ri) acts 
transitively on this orbit. 

Proof: The orbit of any pseudo-pure state p G Dy iV) under U(2n) is homeomorphic 
to CP^"~^ by corollary n The assertion that Sp(n) acts transitively on the unitary 
orbits of pseudo-pure states hence follows directly from the well-known fact that 
Sp(n) acts transitively on CP^"^^ via the isomorphism : H" ^ C'^" discussed 
in [Appendix A| Since the orbit oi p = 5^/2n under U(2n) consists of a single point, 
the assertion that Sp(n) acts transitively on this orbit is trivial. ■ 



Proposition 4 Let p — diag(Ai, • • • , \2n) where Q < \i < 1 are the eigenvalues of p 
counted with multiplicity, satisfying X)i=i ^ ^- Then the orbit of p is a homogeneous 
manifold of real dimension at most 2n^ . 

Proof: Since p is diagonal the isotropy subgroup Gp contains the maximal 
torus T" of Sp(n), namely all matrices of the form diag(zi, • • • , z„, zi, • • • , z„) with 
Zi S U(l) ~ S^. Since Gp is a closed subgroup of the Lie group Sp(n), Proposition 
12 implies that the orbit is a homogeneous manifold of real dimension no more than 
dimR(Sp(n)/r") = n{2n -f 1) - n = 2n^. ■ 

For certain special cases we can improve this bound on the orbit dimensions. 

Proposition 5 If p = diag(cr„, cr„) where (t„ is a diagonal n x n density matrix but 
not a multiple of In then the sympletic orbit of p is a homogeneous manifold of real 
dimension at most 2n^ — 1. If p = diag(a/„, bin) with < a, 6 < 1 and a ^ b then the 
orbit of p under Sp(n) is a homogeneous manifold of real dimension r? -|- n. 

Proof: In the former case, observe that the isotropy subgroup at p contains not 
only the maximal torus T", but all matrices of the form {zJ | z e S*^}. Since the 
matrix J is not symplectic-equivalent to any element of the maximal torus T", the 
first statement follows immediately from the previous proposition. 

In the second case, note that any element Q that belongs to the isotropy subgroup 
at p must satisfy QpQ^ = p. Since Q must also be in Sp(n), it must be of the form 
p5ll for some A, B G endoc(C"'). This implies the matrices A and B must satisfy the 
relations: aAA^-i-bBB'^ = ain and aB-B'^ + MA^ = 6/„ a.nd hence ab{a'^ -b'^)BB'< = 0. 
As < a, 6, < 1 and a 7^ 6 by assumption, we must have B = and A € U(n). Thus, 

the isotropy subgroup at p is ^ A* ^ with A E U(ri). By Theorem the orbit 

of p is a homogeneous manifold of real dimension n{2n + 1) — — n"^ + n. ■ 
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N = 


= 4 


U(4) 


Sp(2) 


P = 


diag(a, a, a, a) 








P = 


diag(a, 5, b, b) 


6 


6 


P = 


diag(a, a, b, b) 


8 


6 


P = 


diag(a, b, c, c) 


10 


< 8 


P = 


diag(a, b, c, d) 


12 


< 8 



N = 


= 6 


U(6) 


Sp(3) 


P = 


diag(a, a, a, a, a, a) 












P = 


diag(a, 6, 6, 5, b, b) 


10 


10 






P = 


diag(a, a, 6, 6, 6, 6) 


16 


< 21 - 


10 


= 11 


P = 


diag(a, 5, c, c, c, c) 


18 


< 21 - 


10 


= 11 


P = 


diag(a, 5, 6, c, c, c) 


22 


< 2 X 32 = 


= 18 


P = 


diag(a, 6, c, d, d, d) 


24 


< 18 






P = 


diag(a, a, b, b, c, c) 


24 


< 18 






P = 


diag(a, a, a, b, b, b) 


18 


<32 + 


3 = 


= 12 


P = 


diag(a, b, c, c, d, d) 


26 


< 18 






P = 


diag(a, 6, c, d, e, e) 


28 


< 18 






P = 


diag(a, 6, c, d, e, /) 


30 


< 18 







Table 2. Dimensions of the orbits of various types of ensembles classified by their 
spectrum under the unitary group U(2n) and the symplectic group Sp{n) for n = 2 
and n = 3. The table shows that the orbits under the symplectic group are in 
general much smaller than the orbits under the unitary group, except in the case 
of pure-state-like (and completely random) ensembles, for which the sympletic 
orbits have the same dimension as the unitary orbit. Note that in all cases it is 
assumed that different letters a,b, . . . represent different values in [0, f ] . 



Proposition 6 If p = diag{D2n-2e, cuhe) , where D2n-2i is a diagonal {2n~2£)x{2n— 
2£) complex matrix, l2e is the identity 21 x 2£ complex matrix, and a G [0, 1], then the 
orbit of p is a homogeneous manifold of real dimension at most n{2n+ 1) — £{2£ + 1). 

Proof: We note the isotropy subgroup at p contains the group 

l2n-2e \ I Q g Sp(i')}, which is isomorphic to the subgroup Sp(^). Hence, 



Q2i 

the orbit of p is a homogeneous manifold of real dimension at most n{2n+l) — £{2£+l) 
as desired. ■ 

In table 121 we compare the dimensions of the orbits of various types of ensembles 
under U(2n) with those of the orbits under Sp(n) for n = 2 and n = 3. The 
classification of the orbits is based on the spectrum of p. Note that the spectrum of p 
uniquely determines the unitary orbit (or equivalence class) p belongs to. However, two 
density matrices with the same spectrum may belong to different symplectic orbits. 
For instance, po = diag(a, 6, a, b) and pi — diag(a, a, 6, b) are unitarily equivalent but 
belong to different orbits under the symplectic group as defined above. See example 
1 in [TT]- 
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5. Conclusion 

We have shown that the action of a Lie group on the set of quantum states endows 
it with the structure of a stratified set with strata given by the orbit manifolds. In 
particular, we studied the stratification of the set of states induced by the action of 
the unitary group, which is especially useful since the unitary orbits are of interest in 
quantum control and computing, where they determine the maximal set of quantum 
states that are reachable from a given set via a coherent control or by applying a 
unitary gate. Furthermore, there are many properties of quantum states such as von- 
Neumann or Renyi entropy that depend only on the unitary orbit the state belongs 
to. It therefore makes sense to define these functions on the unitary stratification. 
We have shown that the unitary orbits can be identified with flag manifolds whose 
type and dimension depend only on the multiplicity of the eigenvalues of the states 
belonging to the orbit. We have also determined the dimensions of the orbit manifolds 
and shown that we can define a partial ordering related to the degree of disorder in 
the system on this stratification via majorization. 

To better understand the geometry and structure of the set of quantum states, we 
studied the embedding of the quantum states and their associated orbit manifolds of an 
n-level system into real Euclidean space provided by the coherence vector. We showed 
that the coherence vector we defined always maps the quantum states into a closed 
ball in M" in such a manner that the orbit manifolds are mapped to submanifolds of 
spheres of fixed radius from the center. For n = 2 this embedding is surjective, hence 
justifying the identification of the set of quantum states with the closed ball in R^, and 
the identification of the orbits with concentric spheres inside this ball. By comparing 
the dimensions of the orbits we also showed that this embedding is no longer surjective 
for n > 2 and the orbit manifolds in this case are proper submanifolds of spheres in 
R" ~^ of lower dimension. The manifold of pure states is always a submanifold of 
the boundary sphere, which contains no other orbits, while all other spheres of fixed 
distance from the center generally contain infinitely many disjoint orbits of varying 
dimensions, and depending on the distance of the sphere from the center, a positive 
measure set of points which do not correspond to quantum states at all. A detailed 
analysis for the three-level case was provided. 

Finally, we studied systems whose natural evolution is restricted to a subgroup 
of the unitary group such as the symplectic group due to dynamical symmetries. 
We showed that we can define a refined stratification based on the smaller orbits of 
this subgroup. In case of the symplectic group we have shown that the orbits of all 
pseudo-pure states agree with the unitary orbits, while the symplectic orbits of all 
other (mixed) states have lower dimension than the unitary orbits. From the point of 
view of control this means that we can control pure and pseudo pure states for such 
systems but not generic ensembles. 
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Appendix A. The symplectic group 

The symplectic group Sp(n) is usually defined as the Lie group of automorphisms on 
H", where H is the skew- field of quaternions, that preserve the canonical symplectic 
inner product 

n 

(q|q')=Eqiqi' 

i=l 

where q and q' are n- vectors whose entries qi are quaternions and conjugation q is 
over H. 

The (skew-field) of quaternions H can be regarded as a vector space over M with 
the standard basis {1,81,62,63} subject to the multiplicative relations: e? = —1 for 
1 < i < 3 and eiej = — ejCi — e^ for any even permutation k) of the set (1, 2, 3). 
Since the field of complex numbers C is isomorphic to M - l©IR-ei and every quaternion 
q can be written as 

q = qo + qiei + q2e2 + qses = (qo + qiei) -|- e2(q2 - qaei), 

we may also regard EI as a vector space over C with basis {1,62}. We therefore 
obtain an isomorphism <j) : H" C^" of complex vector spaces via 4>{qi, •••,(?«) — 
(zi, • • • , Z2n) where qi = Zi + Zn+ie2- Consequently, 

n / 2n \ / n 

(q I q') = I^qiq^ - J2 + ^'M H^^'^n+i - 

i=l \i=l / \i=l 

The isomorphism (j) allows us to identify an isometry A of H", with a complex 
automorphism of C'^" that preserves both the canonical Hermitian inner product and 
the canonical skew-symmetric bilinear form on C^" defined by 

n 

i=l 

where z = (zi, • • • , Z2„) and z' = (z'^, • • • , z'^J. 

Since A preserves the canonical Hermitian inner product on C^", A G U(2n). 
Since A leaves invariant the canonical skew-symmetric bilinear form 5 (•,•), it is equiv- 
alent to having A^JA = J . 
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